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Abstract 

To obtain the most accurate predictions for the Higgs masses in the minimal su- 
persymmetric model (MSSM), one should compute the full set of one-loop radiative 
corrections, resum the large logarithms to all orders, and add the dominant two-loop 
effects. A complete computation following this procedure yields a complex set of for- 
mulae which must be analyzed numerically. We discuss a very simple approximation 
scheme which includes the most important terms from each of the three components 
mentioned above. We estimate that the Higgs masses computed using our scheme 
lie within 2 GeV of their theoretically predicted values over a very large fraction of 
MSSM parameter space. 
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1 Introduction 



With the LEP-2 Collider now beginning operation, a detailed assessment of its discovery 
capabilities has recently been completed The search for the Higgs boson will play a central role 
in the LEP-2 program. A common rule of thumb states that for a center-of-mass energy y/s, a Higgs 
boson with m h o < \fs — 100 GeV can be discovered (at the 5-a level), assuming a sufficient integrated 
luminosity. (Slightly higher masses can be ruled out at 95% CL if no Higgs signal is seen.) This means 
that with yfs = 192 GeV, one can expect to rule out or discover the Higgs boson if its mass is <^ mz- 

In the Standard Model (SM), this Higgs mass discovery reach is not particularly impressive. 
Without any assumption about physics at higher energy scales (above, say, 1 TeV), Higgs masses up 
to about 700 GeV are possible (within the context of a "weakly-coupled" scalar sector). On the other 
hand, if one assumes that no new physics (beyond the SM) enters up to the Planck scale, then one can 
deduce an upper bound of m^o < 175 GeV g. In this light, the LEP search is perhaps more significant. 
But, in the SM, there also exists a lower bound from stability arguments ||. Under the assumption of 
no new physics up to the Planck scale, for mt = 175 GeV, it follows that m h o > 125 GeV. Moreover, 
a Higgs boson discovery at LEP would imply that new physics must enter at an energy scale below 
100 TeV. One of the leading candidates for such new physics is low-energy supersymmetry. 

In the minimal supersymmetric extension of the Standard Model (MSSM), the Higgs sector ^ 
consists of two Higgs doublets of hypercharge ±1. The scalar spectrum consists of two CP-even scalars, 
hP and H° (with m^o < m^o), a CP-odd scalar A and a charged Higgs boson pair H ± . Unlike the 
SM, the Higgs self-couplings are not independent parameters; but they are related to the electroweak 
gauge couplings. As a result, the Higgs mass spectrum is constrained. At tree-level, all Higgs masses 
and couplings depend on two parameters: m A o and the ratio of Higgs vacuum expectation values, 
tan/3. The charged Higgs squared-mass is given by (m 2 H± )o = m 2 A0 + mfy, while the CP-even neutral 
Higgs squared-masses and corresponding mixing angle a are obtained by diagonalizing the 2x2 mass 
matrix (in the hypercharge basis) 

_ ( m-^o sin 2 /3 + m| cos 2 /? - (m 2 A0 + m|) sin cos [3 \ 
\-(m 2 A0 + m|) sin /3 cos /? m 2 A0 cos 2 (3 + m| sin 2 j3 J ' 

where the subscript indicates tree- level quantities. One can then prove that (m h o)o < mz \ cos2/3|. 
If this tree-level inequality were reliable, then LEP-2 would have the mass reach either to discover 
the Higgs boson or rule out the MSSM! When radiative corrections are taken into account, the Higgs 
mass bound increases |5]||. The dominant contribution to this increase is a term of order Am^ ~ 
3g 2 mtlog(M 2 /m 2 )/(8ir 2 m 2 v ), which arises due to an incomplete cancellation of top quark and top 
squark loops (the cancellation would be complete in the limit of exact supersymmetry). 

The radiative corrections to the Higgs mass have been computed by a number of techniques, 
and using a variety of approximations such as effective potential f|j7|,||]9| and diagrammatic methods 
|5|, ^0| , pj] , p^ , p^1 . Complete one-loop diagrammatic computations of the MSSM Higgs masses have been 
presented by a number of groups O]; the resulting expressions are quite complex, and depend on 
all the parameters of the MSSM. Moreover, as noted above, the largest contribution to the one-loop 
radiative corrections is enhanced by a factor of ra\ and grows logarithmically with the top squark mass. 
Thus, higher order radiative corrections can be non-negligible for a large top squark mass, in which case 
the large logarithms must be resummed. The renormalization group (RG) techniques for resumming 
the leading logarithms has been developed by a number of authors [ p^Jl5|JT^1 . This procedure involves 
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integrating a set of coupled partial differential equations. As a result, the numerical evaluation of the 
one-loop Higgs masses is time consuming and not very suitable for detailed phenomenological analyses. 
The primary goal of this paper is to present a simple algorithm that incorporates the effects of the 
RG-improvement and minimizes the size of the two-loop radiative corrections. We then can apply this 
algorithm to a suitable approximation to the full one-loop corrected Higgs masses. 

We present a successive series of approximations to the one-loop corrected Higgs masses, of 
increasing complexity, each one reflecting finer details of the low-energy supersymmetric spectrum. 
Symbolically, 

m 2 H± = (m 2 „±) + (Am%±) . + (Am 2 



H ± - \in H ± j Q t y-un H ± j ill -r y-un H ± 
M 2 = M 2 + AM 2 LL + AM 2 mix , (1.2) 

where the subscript refers to the tree-level result, the subscript ILL refers to the one-loop lead- 
ing logarithmic approximation to the full one-loop calculation, and the subscript mix refers to the 
contributions arising from qL~QR mixing effects of the third generation squarks. The CP-even Higgs 
mass-squared eigenvalues are then given by 



.2 1 



Mfi + M 2 2 ± J [M 2 u ~ M 2 2 } 2 + A{M 2 - ^ 2 



m H°,h° - 

and the corresponding mixing angle, a, is obtained from 

sin 2a = — ; 



l 12) 



;i.3) 



[M 2 u - M 2 2 } 2 + 4(M 2 2 ) 2 

M 2 i - Ml, , . 

cos 2a = ii 22 ( 14 ) 



[Mli ~ M 2 22 \ 2 + 4(M 2 2 y 

In the first (and simplest) approximation, squark mixing effects are neglected and the super- 
symmetric spectrum is characterized by one scale, called Mstjsy- We assume that Mstjsy is sufficiently 
large compared to mz such that logarithmically enhanced terms at one-loop dominate over the non- 
logarithmic terms.0 In this case, the full one-loop corrections (e.g., obtained by a diagrammatic 
technique) are well approximated by the one-loop leading logarithmic approximation. In this approx- 
imation, we neglect all squark mixing effects. Explicit formulae for AA4f LL and (Am^ ± )iLL can be 
found in Appendix A. The second approximation incorporates squark mixing effects. These are likely 
to be significant only in the third generation squark sector. This approximation is parameterized by 
four supersymmetric parameters: MgusY ( a common supersymmetric particle mass) and the third 
generation squark mixing parameters: At, Af, and //. The corresponding formulae for AA 7 ^^ and 
(Am 2 H± ) m ix can be found in Appendix B. 

The dominant contribution to the Higgs mass radiative corrections enters through the exchange 
of the third generation squarks. Thus, our third approximation treats this sector more precisely by 
accounting for non-degenerate top and bottom squark masses. This approximation is characterized 
by seven supersymmetric parameters — the three squark mixing parameters mentioned above, three 
soft-supersymmetry-breaking diagonal squark mass parameters, Mq, Mjj, and Md, and a common 

"If this condition does not hold, then the radiative corrections would constitute only a minor perturbation on the 
tree-level predictions. 
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supersymmetry mass parameter Mgusy which characterizes the masses of the first two generations 
of squarks, the sleptons, the charginos, and the neutralinos. A more precise set of formulae for 
AM-1 LL and AA^^ ix incorporating the detailed squark and slepton mass spectrum, can be found 
in Appendix C. Note that setting Mq = My = Mp = Mstjsy reduces this approximation to the 
previous one. Finally, our fourth approximation incorporates a non-trivial neutralino and chargino 
spectrum. This introduces two additional parameters, Mi and M2 which characterize the color singlet 
gaugino masses. The higgsino masses are determined by (jl, which already enters the analysis through 
the third generation squark mixing effects. The relevant correction terms to AAif LL incorporating a 
non-universal chargino and neutralino spectrum is given in Appendix D. 

Given an approximation to the one-loop Higgs mass as described above we then develop a 
simple algorithm to incorporate the leading effects of RG-improvement. In Section 2, we demonstrate 
that the dominant higher order corrections can be absorbed in the expressions for the one-loop cor- 
rected Higgs masses by a suitable re-definition of m t (and m^). Comparisons between the results of our 
analytic approximation and the results of the numerically integrated renormalization group equations 
(RGEs) demonstrate the domain of validity of our approximations. In Section 3, we include the effects 
of squark mixing. We are able to modify our analytic approximation in a simple way to incorporate 
the main squark mixing effects. Section 4 summarizes the results of this paper. In order to make the 
paper self contained, we collect in a series of appendices the necessary formulae required to implement 
our algorithm. 

This paper was motivated in part by the Higgs Bosons working group of the 1995 LEP-2 
Workshop [0] , which examined in detail the phenomenology of the light CP-even and CP-odd Higgs 
bosons (in the MSSM). The radiative corrections to the light CP-even Higgs mass play a crucial role 
in determining the fraction of MSSM Higgs parameter space accessible to LEP-2 at its maximum 
energy of y/s ~ 192 GeV. As a result, the graphs we present in this paper focus on the light CP-even 
Higgs mass, although our formulae can also be used to compute the radiatively corrected masses of 
the heavy CP-even and charged Higgs bosons. While we were completing this work, Carena and 
collaborators published two papers in which analytic approximations to the radiatively corrected 
MSSM Higgs masses are also developed ||. Their methods and emphasis are somewhat different from 
ours. Nevertheless, the final results are quite similar, and our numerical work (in cases where we have 
compared) typically agree to within 1 GeV in the evaluation of Higgs masses. Some of the results of 
our work have been previously reported in Refs. and p8| , and our formulae have been employed 
in the review of Gunion, Stange, and Willenbrock in Ref. |19| . 

2 RG-improved Higgs Masses — No Squark Mixing 

In this section, we examine the simplest case in which the supersymmetric spectrum is char- 
acterized by a single scale called Mstjsy- We begin with the one- loop leading logarithmic expressions 
for the Higgs squared- masses given in Appendix A. Note that for Mstjsy 3> Tnz, the logarithmically 
enhanced terms appearing in the formulae of Appendix A can potentially spoil the perturbative ex- 
pansion. In this case, it is necessary to perform a RG-improvement which resums the leading logs 
to all orders in perturbation theory. The resulting RG-improved perturbative expansion is better 
behaved and more reliable. The numerical effects of the RG-improvement can be significant for values 
of Mstjsy as low as 500 GeV. 
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Numerical integration of the coupled RGEs is a straightforward but time-consuming process. 
Here, we develop a simple analytic formula that closely reproduces the result of the full numerical 
computation over the parameter regime of interest. Specifically, we are interested in values of Mstjsy 
that lie between 200 GeV and 2 TeV. For values of Mstjsy approaching mz, the leading logarithmic 
corrections are of the same size as non-leading corrections not included in the Appendix A formulae; 
however both are small corrections to the treedevel predictions. Moreover, in this regime, the effects 
of the RG-improvement are insignificant and can be neglected. For values of Mstjsy above 2 TeV, the 
supersymmetry breaking scale is becoming unnaturally large (compared to the scale of electroweak 
symmetry breaking). In this section, we shall apply our technique to the onedoop leading logarithmic 
CP-even Higgs squared-mass matrix, A4 2 LL = Mq + AA4 2 LL . 

The matrix .Mijt depends explicitly on the top quark mass. But, which top-quark mass 
should one use? In a diagrammatic analysis, working in an on-shell scheme, one would use the pole 
mass. The analysis based on RG-running would naturally use the running mass, mt(mt). The choice 
between the pole mass and mt(mt) cannot be decided based on one-loop considerations alone. Since 
the dependence on m>t enters only at one-loop, the distinction between various definitions of m% is a 
two-loop effect. We will return to this distinction later in this section. Dependence on mj will be 
considered below as well, although the numerical distinction among different m b choices is small. 

We can now state a simple analytic formula that incorporates the dominant effects of the 
RG-improvement : 

•MlRG ^ A4 2 ill = Mj LL (m t (fJ,t),m b (fjL b )) , IM = VmMsuSY , M6 = V m z M S vsY ■ (2.5) 

That is, we assert that the numerically integrated RG-improved CP-even Higgs squared-mass matrix, 
•Mf RG , is well approximated by replacing all occurrences of mt and in M. 2 LL (m t , mj) by the 
corresponding running masses.^ Before justifying this assertion, we need formulae for ra^) and 
mt{fj). First, consider m^o = 0(mz)- In this case, at mass scales below Mstjsy j the effective theory 
of the Higgs sector is that of a non-super symmetric two-Higgs-doublet model (2HDM). In this model, 
the quark mass is the product of the Higgs-quark Yukawa coupling (h q ) and the appropriate Higgs 
vacuum expectation value: 

m 6 (/x) = -j=hb{n)vx{n), 

mtip) = —^h t (n)v 2 (ii), (2.6) 

where we employ the normalization vf + v 2 = = (246 GeV) 2 . At scales fi < Mstjsy > we 

employ the one-loop non-supersymmetric RGEs of the 2HDM (see e.g., Ref. |p0|| ) for h b , h t and the 
vacuum expectation values. This yields 



d 2 1 



din// 2 " 16 64 tt 2 

d 2 1 

■ m+ 



Qh 2 b + 2h 2 t - 32g 2 s + p 2 
6h 2 + 2h 2 - 32g 2 - ^g' 2 



m 2 , 



din/, 2 -" 64 vr 2 ' ' ' " ~' " ™* ' (2 ' 7) 

For m^o = O (Mstjsy)) the effective theory of the Higgs sector at mass scales below Mstjsy is that 
of the one-Higgs doublet Standard Model. In this case, we define m q (n) = /ig M (/t)v(/i)/\/2, where 

6 In this section, an overline above a quantity will indicate that the replacement of mt [mt] by mt(fit) [mbifJ'b)] nas 
been made. 
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v(mz) — 246 GeV is the one-Higgs-doublet Standard Model vacuum expectation value. In this case 
eq. ( |2.7D is modified by replacing 6/i 2 + 2/i 2 with 6(/if M ) 2 — 6(hf M ) 2 in the RGE for mf (and interchange 
b and t to obtain the RGE for m 2 ). 

To solve these equations, we also need the evolution equations of g s , and g' . But, an approx- 
imate solution is sufficient for our purposes. Since g' is small, we drop it. We do not neglect the h b 
dependence which may be significant if tan j3 is large. Then, we can iteratively solve eq. ( |2.7| ) to one 
loop by ignoring the /i dependence of the right hand side. We find 

fl-±(a s -i(a b + 3a t )) ln(/i 2 /m 2 ) , m A o ~ 0(m z ) , 
m lW = m lW 4(afH _ ^ M)) ta ^ ^ ^ o(Msusy) (2.8) 

where at = hf/Ait, etc., and all coupling on the right hand side are evaluated at m^. Similarly, 

1 - ~ (a. - ^ (at + 3a b )) ln(^ 2 /m|) , m A o ~ 0{m z ) , 
m b (fi) = m b (m z ) x < ^ .) _ , HM ' n , , 9N (2.9) 



i (a s - A (a SM _ a SM } ) ln (/i 2 /m 2 } f mA0 „ 0(M| 



SUSY J 



For intermediate values of m^o, one may extrapolate the above formulae between the two regions 
indicated in the step function approximation. Using the results of eqs. (2.8) and ( ^li| ) in eq. (2.5), 



and diagonalizing the squared-mass matrix of eq. (2.5) yields our approximation to the RG-improved 
one-loop neutral CP-even Higgs squared-masses. 

Before justifying the above results, we exhibit a numerical comparison among various compu- 
tations of the one- loop corrected light CP-even Higgs mass. First, we evaluate two expressions for the 
RG-unimproved one- loop Higgs mass — the one-loop leading log Higgs mass calculated from A^ 2 LL and 
from a simplified version of -M 2 LL i n which only the dominant terms proportional to mf are kept. In 
the latter case, we denote the CP-even Higgs squared-mass matrix by .M 2 LT = Mq + A7W 2 LT , where 



In many analyses of AA4f LT and AAif LL that have appeared previously in the literature, the 
Higgs mass radiative corrections were evaluated with the pole mass, rat- Some have argued that one 
should take wit to be the running mass evaluated at nit, although to one- loop accuracy, the two choices 
cannot be distinguished. Nevertheless, because the leading radiative effect is proportional to mf, the 
choice of nit in the one- loop formulae is numerically significant, and can lead to differences as large as 
10 GeV in the computed Higgs mass. At the end of this section, we justify the choice of using mt{mt) 
as opposed to na^ ^ (prior to RG-improvement) by invoking information from a two-loop analysis. 
Thus, our numerical results for the light CP-even Higgs mass before RG-improvement is significantly 
lower (when Msusy is large) as compared to the original computations given in the literature, for 
fixed m]? ole . In this paper, we have taken mt(mt) = 166.5 GeV in all of our numerical work. We then 
apply our algorithm for RG-improvement by replacing mt and m b by their running masses evaluated 
at nt and A*6, respectively, as specified in eq. d2.5|) . 

We now show examples for m^o = 1 TeV and two choices of tan j3 in Figs. |l| and ^, and for 
wijsfi = 100 GeV and tan/3 = 20 in Fig. |3[[| Each plot displays five predictions for m h o based on the 

c For m A o — 100 GeV and tan/3 = 1.5, the resulting light Higgs mass lies below experimental Higgs mass bounds 
obtained by the LEP collaborations ]2l|. 
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following methods for computing the Higgs squared-mass matrix: (i) -Mf LT ; (ii) -Mf LL ; (hi) M. 2 ilt (i y ) 
-M 2 ill; and (v) -M 2 RG [see footnote below eq. (^5|) ]. The following general features are noteworthy. 
First, we observe that over the region of MgusY shown, M.\yig — .M 2 ill- in fact, m h o computed from 
-M 2 ill is within 1 GeV of the numerical RG- improved m^o in all sensible regions of the parameter 
space (1 < tan/3 < mt/mb and mt, m^o < Msusy < 2 TeV). For values of Msusy > 2 TeV, the Higgs 
masses obtained from .M 2 ill begin to deviate from the numerically integrated RG-improved result. 
Second, the difference between m h o computed from .Mf LTj and from M.\ RG is non-negligible for large 
values of Mgusy ; neglecting RG- improvement can lead to an overestimate of m^o which in some areas 
of parameter space can be as much as 10 GeV. Finally, note that although the simplest approximation 
of m^o based on .M 2 LT reflects the dominant radiative corrections, it yields the largest overestimate 
of the light Higgs boson mass. 




Fig. 1. The radiatively corrected light CP-even Higgs mass is plotted as a function of Msusy 
for tan f3 = 20 and m^o = 1 TeV. The one-loop leading logarithmic computation [dashed line] 
is compared with the RG-improved result which was obtained by numerical analysis [solid line] 
and by using the simple analytic result given in eq. ( [2.5D [dot-dashed line]. For comparison, 



the results obtained using the leading m\ approximation of eq. ( 2.10 ) [higher dotted line], and 
its RG-improvement [lower dotted line] are also exhibited. Msusy characterizes the scale of 
supersymmetry breaking and can be regarded (approximately) as a common supersymmetric 
scalar mass; squark mixing effects are set to zero. The running top quark mass used in our 
numerical computations is m t (mt) = 166.5 GeV. 



We now give a justification for the prescription proposed in eq. ( |2.5| ) for the radiatively cor- 
rected CP-even Higgs squared-mass matrix. Because the one-loop leading logarithmic correction 
proportional to m\ is dominant, it is sufficient to consider the two- loop leading and next-to-leading 
log contributions proportional to m\a s and mf at- To simplify the presentation we will exhibit the 
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Fig. 2. The radiatively corrected light CP-even Higgs mass is plotted as a function of Mstjsy 
for tan/3 = 1.5 and rriA« = 1 TeV. See the caption to Fig. [|. 
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Fig. 3. The radiatively corrected light CP-even Higgs mass is plotted as a function of Msusy 
for tan/3 = 20 and m^o = 100 GeV. See the caption to Fig. |l|. 
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calculational steps for the case m^o ~ 0(Msusy) 5 n^b = 0, and f3 = n/2. In this case, the effective 
Higgs sector at scales below Msusy is the one-Higgs-doublet Standard Model, which is characterized 
by one Higgs self-coupling parameter A. Supersymmetry fixes the value of A above the supersymmetry 
breaking scale, Msusy, to be 



A(Msusy) = \{g 2 + g' 2 ) (Msusy) 



(2.11) 



The physical Higgs mass is determined at the scale mz (assumed to be much less than Msusy): 
m h o = \(mz)v 2 (mz)- Note that if supersymmetry were unbroken, then A would be fixed to the value 
given in eq. ( [2.11|) , which would imply that m h o = mz- This is the expected result for (3 = it/2 and 
rrijifi > mz [see eq. (l-l)]- In the case of broken supersymmetry, eq. ( 2.11 ) is taken as the boundary 
condition for the one-Higgs doublet Standard Model evolution of A from Msusy down to mz- We 
work in the approximation hf, = g = g' = (a more precise computation is not required here). In this 
approximation, we employ the following one-loop /3-functions 

d . 3 
8^ 



0x = 

Ph 2 = 



din fi 2 
d 

dln/i 



A 



A - h\ - 2Xj v , 



1 



16vr 2 



9 



hi 



8gt 



h 2 , 



where 



1 



v 2 d In /J, 2 



16vr 2 



(2.12) 



(2.13) 



Solving iteratively the RGE for A, and inserting the result in m 2 ^ = X(mz)v 2 (mz), we end up witf 



m 2 h0 



, 3g 2 mj(m t ) 
mz+ M 2 m 2 w ln 



^SUSY 

m} 



1 + 7« + 



Phi 
hi 



hi 



^SUSY 



(2.14) 



where h t = h t {m t ) and m t (m t ) = h t v(m t )/\^2. Clearly, the term in brackets can be absorbed into the 
factor of m\ , thereby converting mt(mt) into the running top quark mass evaluated at a new scale fit- 
Comparison with eq. 



shows that eq. ( 2.14 ) can be rewritten as 
3g 2 mi(fx t ) ( M| USY \ 



m 2 h0 



m 2 z + 



8vr 2 



ln 



rn 



w 



I 5 



(2.15) 



where m = Msusy- Note that in the one-loop leading log analysis presented above, there is 
no distinction between mt{mt) [m h o{mz)\ and the top quark [Higgs boson] pole mass. Moreover, it 
is unclear at what scale to evaluate the factor of mt appearing in the argument of the logarithm in 
eq. ( p. 15 ). These ambiguities can be resolved by a two-loop analysis. 

Ref. [12 1 used a diagrammatic technique to compute the leading two-loop contributions to the 
Higgs pole mass in terms of the top quark pole mass. In the approximation that hf, = g = g' = and 
Msusy 3> irit (and in the absence of squark mixing), the result of this computation can be written in 
the following formF] 



m 2 h0 



m 2 z + 



8vr 2 



m 



(ml 



hi 



-^susy 



^SUSY 
TO? 



4a., 



7T 



cut 



13 _ 

7T 



Note that in the approximation where hb = g = g' 
not run for n < mt- 

e In our two-loop analysis, we have defined g 2 /m^y 
constant fed] and mw is the physical W ± mass. 



1+ V v + t) ln 

(2.16) 

= and A is neglected compared with h t , both A(^t) and v(fi) do 
= i^Gf,, where G M = 1.16639 x 10~ 5 GeV~ 2 is the muon decay 
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In the same approximation as above, the top quark pole mass can be expressed (at one loop) in terms 
of the running top quark mass evaluated at rat 



m 



pole 



m t (m t ) 



1 + 



1,11 

4 



(2.17) 



where mt{mt) is defined below eq. ( 2.14| ).P] Inserting the above result into eq. ( 2.16 ) yields 



m 2 h0 



m| + 



3 5 2 



8vr 2 



m 



■ raj (rat) In 



ir 



-^sus y 
mj 



) 



1 + 



ln / -^SUSY 

I m? 



(2.18) 

Given the Higgs mass computed to two-loop next-to-leading logarithmic accuracy, consider the error 
one makes by using eq. ( |2.14| ), where the two- loop next-to-leading logarithms are neglected. The size 
of the error depends on whether one uses the top quark pole mass or mt(mj). By comparing eqs. ( p.!6| ) 
and ( 2.18| ), we see that the two-loop next-to-leading log term is a 10% correction if the top quark pole 
mass is used, whereas it is a 3% correction if mt(mt) is used.| This is our justification for using the 
running top-quark mass in our one-loop analysis. 

One can also derive eq. ( p. 18 ) using the RG-techniques employed above. Only three modifica- 
tions of our one-loop analysis are required. First, we must distinguish between the Higgs pole mass 
(denoted henceforth by ra h o with no argument) and the running Higgs mass evaluated at m h o. Using 
the results of Sirlin and Zucchini (2f| (in the limit of hf, = g = g' = and A <C /it), 



m 2 h0 



4m^A(m /l o) 



1 f at 

1+- ( — 

7T 



Second, we need only the ht and g s dependent parts of the two loop contribution to (3\. 
of eq. ( 2.12j ) is modified as follows^] [^f| 



(2.19) 
That is, (3\ 



/3a + 



30^-32^3, 



(16vr 2 ) 2 

Third, we must take into account a finite correction between v 2 (m t ) and v 2 = Amfy/g 2 , 



v 2 (m t ) 



Am 2 w 



(2.20) 



(2.21) 



This result can be deduced, for example, from Ref. [12|. Repeating the procedure leading to eq. ( p. 14 ) 
using the modifications indicated in eqs. ( 2.19 )-( 2.21 ), we easily reproduce eq. ( 2.18 ). 

Finally, we can address the question of which top quark mass should appear in the argument 
of the logarithm. Using eq. (2.8) to convert the m 2 in the argument of the logarithm to m 2 (n t ), we 
find 



m 2 



3 5 2 



8vr 2 



m 



( M 2 
j n / JW SUSY 



w 



a 



7T 



l + TT — - 



16 



7T 



(2.22) 



One can check that the sum of the terms in the brackets deviates from one by less than 



■^We caution the reader that Ref. |2^] defines m t (mt) = 2~' i ^ 4 GJ' l 1 ^ 2 ht(mt), which differs slightly from the definition 
of m t (mt) used in our two- loop analysis. 

9 In Ref. jl^] it is shown that the two- loop non- logarithmic piece is negligible. 

h The numerical impact of the two-loop RGEs for the Higgs mass computation was examined in Ref. Jl^] and shown 
to be a few percent of the one-loop radiative corrections. 
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Since the leading mf term provides the dominant source of the neutral Higgs mass radiative 
correction, it follows that our algorithm of replacing m% with mt{^t) in the one-loop formula for -Mf L L 
successfully reproduces the most important aspects of the RG-improvement while minimizing the 
effects of the non-leading-logarithmic two-loop effects, which are numerically small. Note that in the 
numerical work, we use eq. (|2.8| ) to compute the running top-quark mass at the scale of \i% in terms 
of mt(mt). The latter quantity is expressed in terms of the top quark pole mass using the two-loop 



QCD corrections of Ref. |23| (approximately a 6% effect) and the electroweak corrections of Ref. [24] 
(a —0.4% correction for m t = 175 GeV and m h o ~ mz)^ For m^ ole = 175 GeV, this yields roughly 
int(mt) = 166.5 GeV. 

We may also apply our algorithm to the radiatively corrected charged Higgs mass. One can 
check that the twice-iterated two-loop leading log proportional to m 2 m 2 is correctly reproduced. 
However, in contrast to the one- loop radiatively corrected neutral Higgs mass, there are no one-loop 
leading logarithmic corrections to m 2 H ± that are proportional to mf . Thus, we expect that our charged 
Higgs mass approximation will not be quite as reliable as our neutral Higgs mass approximation. 



3 RG-Improved Higgs Masses — Squark Mixing Effects Included 

In the previous section, we focused on the radiative corrections under the assumption that 
one mass scale, Msusy characterizes the super symmetric masses. This is probably not a very realistic 
assumption. In this section, we focus on the effects arising from the mass splittings and qi^qR mixing in 
the third generation squark sector. The latter can generate additional squared-mass shifts proportional 
to mf and thus can have a profound impact on the radiatively corrected Higgs masses.^ In order to 
perform a RG-improvement of these effects, we follow the approach of Ref. fl~6|| . First, we define our 
notation (we follow the conventions of Ref. f27j). In third family notation, the squark mass eigenstates 
are obtained by diagonalizing the following two 2x2 matrices. The top-squark squared- masses are 
eigenvalues of 

/ Mq + m\ + ti,m 2 z m t Xi 



m t X t M 2 +m 2 + t R m 2 z ) ' ^' 2 ^ 

where X t = A t — /i cot [3, = (| — e-t sin 2 0w) cos 2/3 and t R = e t sin 2 6\y cos 2/3, with e t = +2/3. The 
bottom-squark squared-masses are eigenvalues of 

L Q 



Mq + mj + b L m\ m b X b 



m b X b M 2 D + mf + b R m\ 



(3.24) 



where X b = A b — /xtan/3, b^ = (—5 — e b sin 2 9\y) cos 2/3 and 6r = e b sin 2 9\y cos 2/3, with e b = —1/3. 
Mq, Mjj, Md, At, and A b are soft-supersymmetry-breaking parameters, and \x is the supersymmetric 
Higgs mass parameter. We treat the squark mixing perturbatively, assuming that the off-diagonal 
mixing terms are small compared to the diagonal terms. Two levels of approximation are considered. 
At the first level, we take Mq = Mjj = Mr, = MgusYj where Mstjsy is assumed to be large compared 
"Since we are computing the radiatively corrected Higgs mass for arbitrary m A a and tan/3, it is no longer appropriate 



to simply use eq. (2.17) for the top quark pole mass in terms of mt(mt)- 

3 The effects of mass splittings from other supersymmetric sectors to the radiatively corrected Higgs masses are quite 
small and will be neglected in the numerical analysis of this paper. However, in order to be complete, we present in 
Appendix D analytic expressions for the contributions from mass splittings in the chargino and neutralino sectors. 
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to mz- Thus, the radiatively corrected Higgs mass is determined by m^o, tan/3, Mstjsy > A t , A b , and 
fi. At the second level, we allow Mq, My, and Md to be arbitrary. However, as before, we assume 
that these soft-supersymmetry-breaking masses are large compared to mz- 



At one-loop, the effect of the squark mixing is to introduce the shifts A.M mix and (A 



m. 



For Mq = Mjj = Md = Mstjsy 5 the relevant formulae are given in appendix B. For simplicity, we 
focus on this case in the remainder of this section. (For the case of arbitrary Mq, Mjj and Md, see 
Appendix C.) As in Section 2, we note the dependence of A_M mix (m(, mj) on the top and bottom 
quark masses. If we evaluate mt and m;, at a suitable scale, then it is possible to account for the 
dominant effects of the RG-improvement. Thus, for non-zero squark mixing, eq. ( |2.5[ ) is generalized 
to: 

•Mirg ^ -M 2 ill + A_M 2 mix = M\ hh (m t (n t ),m b (fi b )) + AM 2 mix (m t {pLi),m b {^ b )) , (3.25) 



where \i b and \it are defined in eq. (|2.5|) , and fig = Mstjsy (q = t,b). We have extended the overline 
notation of Section 2; note in particular that the scales at which one evaluates nit and m b are different 
in and AA4 mix , respectively. Intuitively, the squark mixing correction arises from integrating 

out heavy squarks that appear in one-loop corrections to Higgs scalar four-point functions. As a result, 
one should choose /ig to coincide with the mass of the heaviest squark. A more detailed justification 
will be given below. 

Following the discussion of Section 2, we compare the value of m h o computed by different 
procedures. Prior to RG-improvement, we first compute m h o by diagonalizing -M 2 LL + A7W mix . Next, 
we perform RG-improvement as in Ref. by numerically integrating the RGEs for the Higgs self 
couplings starting from supersymmetric boundary conditions, and inserting the results into the diag- 
onalized tree-level mass matrix. Finally, we compare the latter with m^o computed by diagonalizing 



M 2 ILL + AM 2 m i x given by eq. (|3^ ). These comparisons are exhibited in a series of figures. First, we 



plot m/jO vs. Xt/MsusY for Mstjsy = "iA° = ~^ = 1 TeV for two choices of tan (3 in Fig. || [tan = 20] 
and Fig. ||| [tan/3 = 1.5]. Note that Fig. ||| is of particular interest, since it allows one to read off the 
maximal allowed value of m^o for Mstjsy < 1 TeV. This maximal value occurs for \Xt\ — 2.4Mstjsy- 
The reader may worry that this value is too large in light of our perturbative treatment of the squark 
mixing. However, comparisons with exact diagrammatic computations confirm that these results are 
trustworthy at least up to the point where the curves reach their maxima. From a more practical 
point of view, such large values of the mixing are not very natural; they cause tremendous splitting in 
the top-squark mass eigenstates and are close to the region of parameter space where the SU(2) xU(l) 
breaking minimum of the scalar potential becomes unstable relative to color and/or electromagnetic 



breaking vacua 28]. 



In Figs. § and [| fi = — 1 TeV, i.e., as X t = A t — /icot/3 varies, so does A t . In fact, for 
m^o 3> mz, the dominant one-loop radiative corrections to m? depend only on X t and Mstjsy [see 



eq. (B.4)1, so that for fixed X t , the /i dependence of m h o is quite weak. This is illustrated in Fig. y. 
For values of m A o ~ 0(mz), the fi dependence is slightly more pronounced (although less so for values 
of tan/3 3> 1) as illustrated in Fig. [?]. We also display m h o as a function of Mstjsy for a number of 
different parameter choices in Fig. [8|. In Fig. ||, we exhibit the tan (3 dependence of m^o for two different 
choices of Xf. Again, we notice that our approximate formula [eq. ( |3.25| )], which is depicted by the 
dot-dashed line, does remarkably well, and never differs from the numerically integrated RG-improved 
value (solid line) by more than 1.5 GeV for Mstjsy < 2 TeV and tan/3 > 1. 
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Fig. 4. The radiatively corrected light CP-even Higgs mass is plotted as a function of 
Xt/MsvSY, where Xt = At — /icot/3, for Msusy = rn^o — — /i — 1 TeV and tan/3 = 20. 
See the caption to Fig. |l|. 
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Fig. 5. The radiatively corrected light CP-even Higgs mass is plotted as a function of 
Xt/MgusY, where X t = A t — fxcot/3, for Msusy = mA° = — M = 1 TeV and tan/3 = 1.5. 
See the caption to Fig. |l|. 
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Fig. 6. The radiatively corrected, RG-improved light CP-even Higgs mass is plotted as a 
function of Xt/MsvSY, where Xt = At — [xcotfi, for MsuSY = tia° = 1 TeV and two choices 
of tan/3 = 1.5 and 20. Three values of fj, are plotted in each case: —1 TeV [dashed], [solid] 
and 1 TeV [dotted]. Here, we have assumed that the diagonal squark squared-masses are 
degenerate: Mq = Mjj = Md = -Msusy- 




Fig. 7. The radiatively corrected, RG-improved light CP-even Higgs mass is plotted as a 
function of Xt/MsuSY for Msusy = 1 TeV and m^o = 100 GeV. See the caption to Fig. || 
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Fig. 8. The radiatively corrected, RG-improved light CP-even Higgs mass is plotted as a 
function of -Msusy for X t = 2.4Msusy for three choices of (tan/3, niA°)= (20,1), (1.5,1), and 
(1.5,0.1), where m^o is specified in TeV units. The solid line depicts the numerically integrated 
result, and the dot-dashed line indicates the result obtained from eq. ( [3.25 ). 




Fig. 9. The radiatively corrected, RG-improved light CP-even Higgs mass is plotted as a 
function of tan (3 for Msusy = 1 TeV and m^o = 250 GeV, for two choices of X t = and 
X t = 2.4M SUSY . See the caption to Fig. |. 
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Finally, we will present a brief derivation of eq. ( |3.25| ). As in Section 2, we simplify the analysis 
by taking m^o ~ C(Msusy) 5 mb = 0> and f3 = tt/2. To incorporate the effects of top-squark mixing, 
eq. ( [2.11 ) is modified as follows 



A (Mj 



SUSY J 



+ 



3A 2 



8tt 2 M s 2 usy 



A\ 



12M| USY 



f4 {Mi 



We may now repeat the steps of Section 2 to obtain: 

2,g 2 miA 2 t 



(Amlo) . 

V / mix 



87r 2 m^M s 2 USY 



1 



12M| USY 



i + 2k + 



hi 



SUSY J 



In (Mksv 
\ mi 



(3.26) 



(3.27) 



Note the additional factor of 2 multiplying the logarithmic term as compared with eq. (2.14). Thus, by 
choosing the running top-quark mass at the scale [ii = Msusy, one completely absorbs the logarithmic 
term, 



mix 8vr 2 m 2 y M| USY 



1 



A 2 



12Mg UgY 



(3.28) 



4 Discussion and Conclusions 



The masses and couplings of the Higgs bosons of the minimal supersymmetric model (MSSM) 
are calculable in terms of two parameters (m^o and tan (5) at tree-level. Including radiative corrections 
brings in a dependence on the MSSM particle spectrum, with dominant sensitivity to the top quark 
mass and the top squark masses and mixing parameters. Perhaps the most important prediction of the 
MSSM Higgs sector is the value of the lightest CP-even Higgs mass, m^o, which is predicted to be no 
greater than mz\ cos 2(3\ + Am, where Am depends on the size of the radiative corrections. This implies 
that a significant fraction of the MSSM Higgs parameter space is accessible to LEP-2. In this regard, 
an accurate computation of m h o is crucial for a reliable assessment of the capabilities of LEP-2 (as a 
function of its center-of-mass energy) . Presently known techniques allow one to compute m^o with an 
accuracy of roughly 2 GeV. If the top squark masses are large (> 500 GeV), then it is crucial to include 
renormalization group (RG) improvement to obtain the accuracy just quoted. In this paper, we have 
presented a simple algorithm for including RG-improvement in the computation of Higgs masses. Our 
proposal can be applied to either the exact one-loop diagrammatic computation, or to approximations 
thereof. If we denote the one-loop CP-even Higgs squared-mass matrix by M 2 (mt, m b , Xt,X b ), where 
we have indicated the dependence on the third generation quark mass and squark mixing parameters, 
then, the RG-improved squared-mass matrix is well approximated by 



M 2 = M 2 [m t {n~),m h {iii),X t ,X h ] - M 2 [m t {^),m b {^ b ), 0, 0] + M 2 [m t ^ t ),m b (n b ), 0, 0] (4.29) 



This summarizes the results of eq. (|2.5|) and eq. ( 3.25| ). Note that eq. ( [4.29D can be applied even in 
the case where the effects of squark mixing are not explicitly separated out. As an example, consider 
the full diagrammatic one-loop computation, whose result we shall denote by -A/f 2 LC . To carry out the 
RG-improvement, one would subtract off the one- loop leading logarithmic pieces, -M 2 LL , and replace it 
with its RG-improved value, _M 2 RG . If we implement the RG-improvement using our analytic method 
described by eq. ( 4.29| ), we obtain 



M 2 1LC = M( LC + M 2 lhL + AM 2 , 



M 2 



ILL 



AM' 



(4.30) 
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In cases that we have checked, we find that A4f LC — -Mill — A-^mix produces no more than about 
a 1 GeV shift in the predicted Higgs masses. Thus, we expect that the numerical results obtained in 
this paper will shift by no more than ~ 1 GeV if non-leading logarithmic terms are included. 

We find that if our algorithm is applied to the leading log one-loop corrections (given in 
Appendices A and C), plus the leading terms resulting from squark mixing (given in Appendices B and 
C), we reproduce the full (numerically integrated) RG- improved value of m h o, to within an accuracy 
of less than 2 GeV (assuming that supersymmetric particle masses lie below 2 TeV). Although we have 
not focused our attention to the heavier Higgs states, we assert that our algorithm also yields accurate 
results for the mass of the heavier CP-even Higgs boson, m H o. Our approximation to the radiatively 
corrected charged Higgs mass is slightly less accurate only because the leading mt enhanced terms are 
not as dominant as in the neutral Higgs sector.^] Thus, the methods described in this paper provide 
a simple and powerful technique for obtaining an accurate determination of the radiatively corrected 
MSSM Higgs masses and CP-even mixing angle. 
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Appendix A: One-Loop Leading Logarithmic Corrections to the MSSM Higgs Masses 

We suppose that the supersymmetric particle spectrum is characterized by one mass scale, called 
-^SUSY (assumed to be larger than m t ) Specifically, we assume that terms proportional to 
ln(Mg USY /Mj^), where Mx stands for the mass of any supersymmetric particle, are much smaller 
than terms proportional to ln(Mg USY /m^) and can be neglected]] This assumption is made to keep 
the formulae that follow relatively simple. In fact, one can account for some spread in the super- 
symmetric mass spectrum by interpreting Mgusy appropriately depending on where it occurs in the 
formulae below. In this case additional correction terms arise, which are treated in Appendices C 
and D. 

In the one-loop leading logarithmic approximation, the CP-even Higgs squared-mass matrix is 
given by: 

M 2 LL = M 2 + AM 2 LL , (A.l) 
where Mq is the tree level result given in eq. ( |1.1| ), and 



(AA^fl)lLL 



g 2 m 2 z cos 2 (3 
96-7T 2 cos 2 6w 



I M 2 



mi 



mi 



+ (l2N c . ; 
\ m % cos 4 p 
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m z cos z p 



+ P b + P f +P g + P 2H In 



^SUSY 
m\ 



+ Q(m A o - m z )(P\H ~ P2H) In 



m 



A" 



mi 



(AM 2 22 ) lhL 



g m z sin (3 
967r 2 cos 2 9w 

+ ( 12N r 
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m% 
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m z sm p 



+ l > ) In ( Mksv 



mi 



+ Q(m A o - m z ){P\H ~ P2h) In 



in 



A" 



mi 



+ 2N C 



mi 



m\ sin (5 



(AM 2 12 ) 1LL 



gm 2 z sin (3 cos (3 
967r 2 cos 2 Ow 

+ (p b - m 



Pt - SN C 



mj 



m? z sin 2 (3 



\ mi 



mi 



m z cos 2 (3 



+ P f + P' + P' 



2H 



l n ( MuSX 
\ mi 



Q(m A o - m z )(PiH + P 2 h) ln 



hi 



A 



mi 



(A.2) 



'For a widely spread supersymmetric particle spectrum where this approximation is inappropriate, these corrections 
are listed in Appendices B and C (squarks and sleptons) and D (charginos and neutralinos) . 
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In eq. ( A. 2 ), the constants Pj are denned as 



Pt 

Pb 

Pf 

P 9 
P' 

9 

P2H 

P' 

r 2H 

PlH 



= N c (1 - Ae t sin 2 W + 8e 2 sin 4 9 



w ) 



N c (l + 4e 6 sin 2 9 W + 8e 2 sin 4 9 W ) 



= N c (N g -l) 



2-4 sin 2 9 W + 8(e 2 + e 2 ) sin 4 9 W 
-44 + 106 sin 2 6 W - 62 sin 4 % , 
10 + 34 sin 2 W - 26 sin 4 9 W , 
-10 + 2 sin 2 9 w -2 sin 4 6 W , 
8-22 sin 2 % + 10 sin 4 9 W , 
-9 cos 4 2(3 + (l - 2 sin 2 % + 2 sin 4 6> w ) cos 2 2/3 . 



+ iV 



2-4sin 2 6> w + 8 sin 4 W 



(A.3) 



with the number of colors, N c = 3, the number of generations, N g = 3, and the quark charges given 
by &i = 2/3 and e& = —1/3. The origin of the various terms above can be found in Ref. The 
most important contribution to the one- loop radiative corrections above arises from the term in AA4 22 
that is proportional to m\. If one discards all other one-loop radiative corrections, then one obtains 



A7W 2 LT given in eq. ( |2.1C| ). 

In principle, the two cases m A o < 0(mz) and m A o 3> mz must be treated separately. However, 
by including the terms proportional to ln(m A o / w-z) in eq. ( A.2[ ) ,p| we interpolate between the small 
and large m A o cases. This yields the correct leading log expression for m h o . The large m A o analysis is 
more conveniently done in a basis where h° and H° are approximate eigenstates. If 1Z is the rotation 
matrix which rotates a 2-vector counterclockwise by angle (3, then TZ~ 1 A4 2 TZ yields 



m 2 



m 



hH 



hH 



m 



H° 



M 2 xl cos 2 (3 + M 22 sin 2 (3 + M\ 2 sin 2(3 \ (M 2 22 - M 2 n ) sin 2(3 + M\ 2 cos 2/3 



\{M 2 22 



Mfi) sin 2(3 + M{ 2 cos 2(3 M 2 U sin 2 (3 + M l 22 cos 2 (3 



M\ 2 sin 2(3 



where the identification of the diagonal elements above follows from the fact that m\ H 
the limit of m A o mz- Using the above results, it follows that for m A o » mz, 
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(A.5) 



In order to compare the results of our analytic approximation to the RG-improved Higgs squared-mass 
[eq. (|2.5|) 1, one must integrate the Higgs self-couplings of the two-Higgs-doublet model from MgusY 



m Strictly speaking, these terms are meaningful for mz >C m A o < Msusy- The step function in eq. (A. 2) implies that 
these terms should not be included if m A o < mz- 
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down to m A o, match onto the one-Higgs-doublet model, and finally evolve the Standard Model Higgs 
self-coupling from m A o down to mz- However, we note that the contribution of the terms proportional 
to \n{m 2 A0 / m 2 z ) are quite small, leading to a mass shift in m h o of no more than 1 GeV (for m A o < 
1 TeV). As a result, we simplified the analysis in generating the graphs shown in Section 2 and 3 by 
omitting the m A o dependence in the radiative corrections.^ That is, we omitted all terms proportional 
to ]n(m A0 /m z ) that appear in the formulae above. We then compared the analytic approximation so 
obtained with the two-Higgs-doublet model numerically integrated from Mgusy down to mz, with no 
matching onto the one-Higgs-doublet model as described above. In this approximation, the dependence 
of the Higgs masses on m A o enters only through the tree-level expressions. The impact of the m A o 
dependence of the radiative corrections will be discussed briefly at the end of Appendix C. 

Although the one-loop leading log terms in m 2 are correctly reproduced by using the results 



given m eqs. (Q -dAH, the analo gous terms in m 2 H0 obtained by the above procedure do not 



catch all terms proportional to ln(m^ 



/m|) 



The reason is related to the fact that below the scale 
m A o, the effective low-energy Higgs sector consists of just one scalar, hP. Thus, for example, for 

^2 _ ™2 i ™2 ™2 . 



Msvsy one would find m H0 



m 



A" 



+ m z sin 2/3, where the parameters on the right hand 



side should be considered as running parameters evaluated at the scale m A o (and similarly for m hH = 
—m\ sin 2/3 cos 2/3). Then, the correct leading log pieces in the unspecified 0{g 2 m 2 z ) terms in eq. ( A. 5 ) 



would be obtained once the running parameters were expressed in terms of physical parameters. We 
shall not provide explicit expressions for the resulting 0{g 2 m 2 z ) terms, since these corrections^ yield 
only a small fractional shift to m 2 H0 . One subtle point involves the definition of tan/3. The angle (3 
that appears in the expression for m 2 in eq. ( |A,5| ) is implicitly defined at the scale mz- However, 
for m A o mz, the Higgs sector is effectively the one-doublet scalar sector of the Standard Model 
below the scale m A o. Thus, another choice is to define tan (5 as the ratio of vacuum expectation values 
evaluated at m A o. To convert between definitions is a simple task. In the expression for m 2 ^ given 
in eq. ( |A.5| ), it suffices to replace the tree- level contribution, m? z cos 2 2/3 (mz), using the leading log 
expression |l6| 



cos 2 2j3{mz) = cos 2 2/3(m^<>) + 



g 2 N c cos 2/3 



m 
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m, cos' 



: /31n 



in 
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.4" 



mi 



. (A.6) 



However, we alert the reader that we have not made this replacement in any of our formulae. If not 
explicitly indicated, all quoted formulae in this paper assume that tan/3 is defined at mz- 

We can improve the above formulae by reinterpreting the meaning of MgusY- For example, all 
terms proportional to ln(M| USY /m 2 ) arise from diagrams with loops involving the top quark and top- 
squarks. Explicit diagrammatic computations then show that we can reinterpret Mf USY = M f Mr . 
Note that with this reinterpretation of -^|usy» the top quark and top squark loop contributions 
to the Higgs masses cancel exactly when Mr = Mr = mt, as required in the super symmetric limit. 



Likewise, in terms proportional to Pb or powers of multiplied by ln(M| ugY / 'm 2 ?) '> we ma Y reinterpret 
MgusY = -^Si^b " Terms proportional to Pf ln(M^jj S Y / m z ) come from loops of lighter quarks and 
leptons (and their supersymmetric partners) in an obvious way, and the corresponding Mg USY can 
be reinterpreted accordingly. Additional contributions arising in the case of non-degenerate squark 

"In producing the graphs shown in Sections 2 and 3, we also have omitted the non-leading logarithmic term propor- 
tional to that appears in AM \ 2 , since this term is not picked up by the numerical integration of the RGEs. 

"Note that the above arguments imply that if m A o = Msusy then the leading mf contribution to the (D(g 2 m%) terms 
is absent, while for mt < m A o < Msusy, the leading mf term is proportional to m\ ln(Mg VSY /m A o ). 
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masses are treated in Appendix C. The remaining leading logarithmic terms arise from gauge and 
Higgs boson loops and their supersymmetric partners. The best we can do in the above formulae 
is to interpret MgusY as an average neutralino and chargino mass. To incorporate thresholds more 
precisely requires a more complicated version of eq. ( |A.2| ), which can be easily derived from formulae 



given in ref. [16|. The corresponding expressions are summarized in Appendix D. However, the impact 
of these corrections are no more important than the non-leading logarithmic terms which have been 
discarded. The most significant non-logarithmic correction is one term that scales with m 2 and thus 



has been included in eq. (A. 2). This is surely the largest of such corrections. One can check that it 
yields at most a 1 GeV shift in the computed Higgs masses. 

Finally, we give the one-loop leading logarithmic expression for the charged Higgs mass. First, 
assuming that m A o < 0(mz), 
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(A.7) 



For m A o 3> mz, additional terms arise proportional to \n.{m 2 A0 / 'my/) . For example, for m A o = Msusy, 

+ m^y, where m 2 ^ is treated 



A 



the leading logarithmic corrections can be obtained from m H± - 
as a running parameter evaluated at m A o . Re-expressing m\y (m A o ) in terms of the physical W mass 
yields the correct one-loop leading log corrections for mz < m A o < Msusy- Again, we omit the 
explicit expressions since these corrections generate only a small relative shift to the heavy charged 
Higgs mass. 

A good approximation to the RG-improved Higgs squared-mass corrections is implemented 
according to the algorithm of Section 2 by replacing mj and m& in the above formulae by the corre- 
sponding running parameters evaluated at fit = v^t-^SUSY and fib = \Jmz Msusy > respectively. 



Appendix B: Leading Squark Mixing Corrections to the MSSM Higgs Masses 

When squark mixing effects are taken into account, the CP-even Higgs squared-mass matrix and rn 2 H ± 
receive additional one-loop corrections beyond those given in Appendix A. We need only consider the 
effect of mixing among the third-generation squarks (since all such mixing effects are proportional to 
the corresponding quark mass). In order to keep the formulae simple, we assume that the diagonal 
elements of the 2x2 squark squared-mass matrices [eqs. ( ^23|) and ( gjg) ] are degenerate and given 



by M| USY . (In Appendix C, we treat the more complex case of non-universal squark squared- masses.) 
It is convenient to define 

X t = A t - fi cot P , Y t = A t + fi tan P , 

X b = A b -fitanP, Y b = A b + ficotp . (B.l) 
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We assume that the mixing terms mtX t and m b X b are not too large|] Then, the elements of the 
CP-even Higgs squared-mass matrix are given by: 



M 2 = Mj LL + AM; 



(B.2) 



where -A/ff LL has been given in eqs. ( |OD and ( gjg ), and 
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(B.3) 



If mz <C m^o < MgusY) then it is again convenient to rotate the neutral scalar Higgs basis 
as discussed in Appendix A. One then finds that the Higgs squared-masses obtained in eq. ( |A.5| ) are 
shifted by the following expressions 
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Squark mixing effects also lead to modifications of the charged Higgs squared-mass. One finds 
that the charged Higgs squared-mass obtained in eq. (A. 7) is shifted by 
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^Formally, the expressions given in this Appendix are the results of an expansion in the variable (Mj 2 — Mj 
where Mi , M| are the squared-mass eigenvalues of the squark mass matrix. Thus, we demand that mtX, 
For example, for Msusy = 1 TeV, values of X t /Msusy Si 3 should yield an acceptable approximation 
formulae presented here. 



•■)/(M?+Ml), 
/M| USY « 1. 
based on the 
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A good approximation to the RG-improved Higgs squared-mass corrections is implemented 
according to the algorithm of Section 3 by replacing mt and m\, in the above formulae by the corre- 
sponding running parameters evaluated at fx = MgusY- 



Appendix C: Non-Universal Squark Mass Corrections 



The formulae of Appendices A and B were obtained under the assumption that Mq = M\j = Mp = 
MgusY) where Mq, M\j, and Mp are the diagonal soft-supersymmetry-breaking squark mass param- 
eters defined in eqs. ( gjp and (|]2D. In this Appendix, we generalize the results of the previous two 
appendices to allow for unequal diagonal scalar masses. 

First, the results of Appendix A are modified as follows. For the neutral CP-even squared-mass 
matrix, replace in eqs. ( |A.2[ )-( [A5|) all occurrences of ln(Mg UgY /m 2 ) with \h(MqMjj / m 2 ), and replace 
all occurrences of ln(M| ugY /m 2 z ) which multiply either mf, m\ or by ln(MgM£)/m|). In addition, 



eq. ( |A.2| ) must be modified by adding the following squared-mass shifts: 
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added to (m? )iLL given in eq. (|A. 



In the limit of m^o 3> mz [using eq. ( |A.4| )], we find that the following squared-mass shift must be 
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For m H ±, replace ln(Mg USY /m 2 ) in eq. O with ipnlMgM^/m?) + ln(M Q M D )/m 2 )]. In 
addition, the following squared-mass shift must be added to (m 2 H± )u J L given in eq. (|A.7|) : 
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where the function g(a, b) is defined by: 

g(a,b) = 



(a - bf 



2 _a_ + b (a 
a — b \b 



(C.4) 



Note that g(a,a) = — l/6a 2 , so that in the limit of Mq = M\j = Mb, all mass shifts given above 
vanish. 

Next, we consider the modifications of the results of Appendix B when the diagonal soft- 
supersymmetry-breaking squark mass parameters are non-degenerate. The shifts in the squared-mass 
matrix elements of the CP-even Higgs boson given in eq. ( |B.3[) are replaced by the following expressions: 
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where the functions B, h, p b and pt are defined as follows: 
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g(a,b) is given in eq. ( |C.4| ) and f(a,b) is given by: 
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The contributions to AA^ 2 ^ proportional to or mf were first given in Ref. 
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To check that eq. ( |C.5|) reduces to eq. ( [B.3D in the limit of Mq = Mjj = Mp, one makes use 
of B(a, a) = l/6a, /(a, a) = — l/2a, g(a, a) = — l/6a 2 , and h(a, a) = 1/a. In the limit of m A o S> mz, 
eq. ( [B.4D is replaced by: 
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The shift in the charged Higgs squared-mass given in eq. 



is replaced by: 
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where the new functions that appear depend on three variables: 
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To check that eq. ( p.9| ) reduces to eq. ( |B.5[ ) in the limit of Mq = Mjj = Md, one should note that 
F(a, c, c) = f(c, a) and G(a, c, c) = —g(a, c). To take the limit of a = c, we may use the results quoted 
below eq. ( p.7|) . 

The above results correspond to the leading terms of an expansion in (M 2 — M|)/(M 2 + M 2 ), 
where M 2 , M 2 are the squared-mass eigenvalues of the squark mass matrix. If the squark squared-mass 
splitting is large, one needs a slightly better result. This is easily obtained from the formulae above 
for the CP-even Higgs squared-masses as follows. Terms multiplying factors of m b involve functions 
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of the arguments Mq and Mp. In these terms, replace Mq and M|> by the corresponding bottom 
squark squared-masses. Similarly, terms multiplying factors of mt involve functions of the arguments 
Mq and M^j. In these terms, replace Mq and My by the corresponding top squark squared-masses. 
One can check that this rule correctly reproduces the terms in the CP-even Higgs squared-mass matrix 
proportional to mf and vn\ obtained by the effective potential methods of Ref. ||]. Moreover, we have 
explicitly verified that the terms proportional to are also correctly reproduced at (3 = it/2 by 
comparing with the exact one-loop computation of Ref. ||. In the case of the charged Higgs mass, 
the improvement required in the case of large squark squared-mass splitting is more complicated. In 
particular, the simple rule quoted above does not apply, since there are functions appearing in the 
charged Higgs mass formulae that contain both My and Mp as arguments, so the effects of top and 
bottom squarks do not factorize. 

Finally, we must address the question of effective scales fj, q and fig (q = t, b) introduced in 
Sections 2 and 3, respectively, to account for the renormalization group improvement. By iteratively 
solving the RGEs with two different squark masses we can show that the two-loop leading log correc- 
tions can be absorbed into the one-loop expression by choosing 



j4_ = ln 2 (M h /m t )+ln 2 (M h /m t ) 
m\ " \n{M~ ti M i2 /ml 



and similarly for [with m< replaced by mz in the arguments of the logarithms] . Similarly, one can 
modify the analysis at the end of Section 3 to show that the fig defined below eq. ( 3.25[ ) generalizes 



to Hq = Mq 2 (q = t,b), i.e., the mass of the corresponding heaviest squark mass eigenvalue. 

To illustrate the effects of the corrections described in this section, we exhibit in Fig. 10 the 
radiatively corrected RG-improved light CP-even Higgs mass as function of Msusy for m^o = At = 
—fi = Mq = M\j = Md = Msusy and two choices of tan/3 = 1.5 and tan/? = 20. The graphs have 
been computed using the algorithm of eq. ( |3.25 ). The dotted line is based on the formulae presented 



in Appendix A and B, where terms proportional to ln(m^ / 'rn%) and the non-leading logarithmic term 
proportional to m\ have been omitted, as described below eq. ( |A.5| ). The importance of these two 
terms is illustrated by including them in the dashed curve of Fig. 110. This curve also includes the effect 
of replacing M| USY in the leading logs of Appendix A by the appropriate product of third generation 
squark masses, as described below eq. (A. 6). We call this the "improved one-scale" approximation. 



Finally, we improve on this approximation by employing the correction terms described in this section, 
with the replacement of Mq, Mjj, and Md by the corresponding third generation squark masses as 
described below eq. ( |C.11| ). The result of this improved treatment of third generation squark thresholds 
is depicted by the solid line in Fig. [H]. One can see that the improved treatment of thresholds typically 
introduces no more than a 2 GeV mass shift in the predicted value of m h o. 

It is important to emphasize that all the formulae given in Appendices B and C are based on 
the assumption that the top and bottom squark masses are substantially larger than mz- Fortunately, 
in most cases of interest, the error introduced if the squark masses are not sufficiently heavy is 
rather small, of the same size as the non-logarithmic terms that have been systematically neglected 
in this paper. Thus, one can confidently use the formulae given in this paper over most of the 
relevant supersymmetric parameter space, and expect an accuracy in the computed Higgs masses 
within ~ 2 GeV of their actual values. 
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Fig. 10. The radiatively corrected light RG-improved CP-even Higgs mass is plotted as a 
function of .Msusy for m^o = A t = —fi = Mq = Mjj = Mo = Msusy- Two choices of 
tan/? = 1.5 and tan/3 = 20 are shown. The results are based on the analytic approximation 
given in eq. (|3.25 ). In the one-scale approximation (dotted line), the distinction between 
-Msusy and third generation squark masses are neglected. In addition, radiative corrections 
proportional to \i\.{m 2 A0 / m 2 z ) and the non-leading log term proportional to m 2 are neglected, as 
discussed below eq. (A.5). These approximations were also made in the graphs of all previous 
figures. In the improved one-scale approximation (dashed line) , the latter two neglected terms 
are included. In addition, the factors of -M| USY in the argument of the leading logs of Appendix 
A are replaced by the appropriate product of third generation squark masses, as described below 
eq. ( A.6). In the improved threshold approximation, we include the corrections of Appendix 
C with the replacement of Mq, Mjj, and Mjj by the corresponding third generation squark 
masses as described below eq. (C.ll). 
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Appendix D: Non-universal Gaugino and Higgsino Mass Corrections 



The formulae of Appendix A were obtained under the assumption that the chargino and neutralino 
masses are all degenerate and equal to MgusY (assumed large compared to mz)- In this appendix, 
we relax the assumption of degeneracy of masses, although we still assume that all the chargino 
and neutralino masses are large compared to mz- If the latter condition is not satisfied, then the 
logarithmic pieces that we keep are of the same order of magnitude as the non-logarithmic pieces that 
we omit. Nevertheless, the contributions of the charginos and neutralinos to the Higgs masses are 
small (never exceeding 2 GeV), so our approximations will be accurate at this level over the entire 



neutralino and chargino parameter space. The results below have been extracted from ref. [16]. 

In the limit where all chargino and neutralino masses are large compared to mz, the chargino 
spectrum consists of two states of mass M 2 and \fi\ respectively, while the neutralino spectrum consists 
of two states of mass Mi, M2, and two approximately degenerate states of mass Here, Mi and M2 
are Majorana masses for the uncolored gauginos, and fi is the supersymmetric Higgs mass parameter 
that also appears in the off-diagonal squark mixing matrix.^] It will be convenient to introduce three 
mass parameters which are equal to either Mi, M2, or \fx\ depending on their relative magnitude: 

111 = max{|/x|, Mi} 
/i 2 = max{|//|, M 2 } 

Hl2 =max{ H, Mi,M 2 } (D.l) 

We consider the case of unequal Mi, M 2 , and \fi\, all of which are assumed to be large compared 
to mz- Then, eq. ( |A.2j ) of Appendix A is modified by adding the following squared- mass shifts: 
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where the Pj are defined in eq. ( |A.3| ). In the limit of m^o 3> mz [using eq. ( A. 4 )], we find that the 
following squared- mass shift must be added to (m^ )iLL given in eq. (| 
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Note that the effect of the terms proportional to ln(Mg USY /m 2 ? ) is to simply remove the MgusY 



dependence in eqs. ( |A.2| ) and ( |A.4| ) that arises from the gauge/Higgs/gaugino/higgsino contributions. 
These contributions are now more accurately described by logarithmic factors that are sensitive to the 
gaugino/higgsino spectrum. 



For m H± , the following squared- mass shift must be added to (mi 1 ±)ii,x J given in eq. (|A.7|) : 
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As above, the effect of the terms proportional to \n(M§^ Y / m \y) 1S to simply remove the MgusY 
dependence in eq. ( |A7D that arises from the gauge/Higgs/gaugino/higgsino contributions. 

One can easily check that if Mi = M2 = = MgusY; then all the mass shifts in this appendix 
vanish. In the case of the neutral [charged] Higgs mass computation, if any of the gaugino/higgsino 
mass parameters is less than mz [mw], then one must remove the corresponding logarithmic term from 
the above expressions. Note that the supersymmetric limit corresponds to all three gaugino/higgsino 
mass parameters zero and tan/3 = 1, in which case two charginos are degenerate with the and 
H^, two neutralinos are degenerate with the Z and H°, one neutralino is degenerate with the photon, 
and one neutralino is degenerate with the (massless) h° and A . In this case, we must remove all 
logarithmic terms above, except for the log terms containing Afg USY . The effect of the latter is to 
precisely cancel the contributions of the gauge and Higgs loops to the CP-even Higgs squared-mass 
matrix [eqs. ( |A.2| ) and (A. 4)] and the charged Higgs mass [eq. (A. 7)]. This cancellation is a consequence 
of unbroken supersymmetry. Of course, as emphasized above, when some of the gaugino/higgsino mass 
parameters are of order mz, then the above corrections terms are small in magnitude and of the same 
order as non- logarithmic corrections not included in this analysis. 
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